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Abstract 

We report on a novel response to biases in deterministic superdiffusion. For its reduced map, we 
show using infinite ergodic theory that the time-averaged velocity (TAV) is intrinsicahy random 
and its distribution obeys the generahzed arc-sine distribution. A distributional limit theorem 
indicates that the TAV response to a bias appears in the distribution, which is an example of what 
we term a distributional response induced by a bias. Although this response in single trajectories 
is intrinsically random, the ensemble-averaged TAV response is linear. 
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Introduction. — Intrinsic randomness in macroscopic observables has been found in a broad 
range ofprocesses more recently from biological transport to fluorescence of single nanocrys- 
tals 1- J Diffusion coefficients in biological transports show large fluctuations [l, 2|. The 
ratio of an on-state in the fluorescence of nanocrystals does not converge to a constant and 
is different in each quantum dot j^. The randomness of time-averaged observables can be 
characterized by power-law trapping time distributions in stochastic models, such as those 
in continuous time random walks (CTRWs) 0, and dichotomous stochastic processes j^. 
Such randomness is due to the breakdown of the law of large numbers stemming from a 
diverging mean trapping time. 

Dynamical systems with infinite-invariant measures can be viewed as stochastic processes 
generating random time-averaged observables. Infinite ergodic theory plays an important 
role in elucidating such observables [lUl. It guarantees that a time average of an obser- 
vation function converges in distribution. In other words, time-averaged observables are 
intrinsically random if the invariant measure cannot be normalized. In dynamical systems 
generating subdiffusion, the distributional limit theorem for the diffusion coefficients ob- 
tained by the time-averaged mean square displacements (TAMSDs) has been shown using 
infinite ergodic theory 12 1. 

In anomalous diffusion, the mean square displacement (MSD) grows non-linearly with 
time, {x{tY) oc (a 7^ 1). Diffusion is called subdiffusion if a < 1 and superdiffusion 
if a > 1. One mechanism generating subdiffusion is the divergence in the mean trapping 
time describing random walks of CTRW. If the mean trapping time diverges, the diffusion 



121]. Three different mechanisms 



coefficients obtained from TAMSDs become random [4, |5|, 
underlying superdiffusion have been identified: One stems from positive correlations in ran- 
dom walks, modeled by a fractional Brownian motion [isl; The second from a persistent 
motions in random walks, called Levy walks [l^; The third from very long jumps in random 
walks, called Levy flights isj. In Levy walks and flights, the second moment of the length 
of persistence motion and of iumps diverges because of a power law. Such power laws are 

n n n 

observed in Hamiltonian systems |16|. rotating flow [17[, polymer diffusion [15|, biological 



transport [18|, intermittent search [19[, and light diffusion (20[. Persistent times in Levy 
walks and trapping times in CTRWs are characterized by indifferent fixed points in deter- 
ministic models |2l|. Therefore, random transport coefficients will be observed in Levy walk 
with the divergent mean persistent time. 



Although random transport coefficients are universal in subdiffusion because of power-law 
trapping times, it is not clear whether time-averaged observables are intrinsically random in 
super diffusions. Moreover, little is known about responses of time-averaged drifts (TADs) 
to biases in superdiffusions whereas in anomalous diffusion a generalized Einstein relation 
holds 22|. In this Letter, we show a distributional limit theorem for TADs using determin- 
istic superdiffusion models related to Levy walks under bias and no bias. Surprisingly, a 
TAD is intrinsically random whether biased and unbiased. The result leads to a distribu- 
tional response to bias, i.e., the response to a bias would be characterized by a change in 
distribution. 

Model. — Anomalous superdiffusion that originates from a persistent motion or a long 
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24| . To study a response to 



jump has been studied for chaotic dynamical systems 
a bias in deterministic superdiffusions related to Levy walks, we propose an asymmetric 
deterministic diffusion model, constructed by introducing an asymmetry in the Geisel model 



23(1 . In particular, we consider the following map T : M — )■ M, 

Xn+l = T{Xn), (1) 

which have a translational symmetry [L = Z) 

T{x + L)= T{x) + L, (2) 

and the map T{x) is given by 

f {x-L) + + L - 1, X e [L, L + c), 

T{x) = (3) 
({x-L)- + L + 1, X G [L - c, L), 

where c is a parameter characterizing an asymmetry. This model corresponds to a Levy walk 

where leftward and rightward persistent time distributions have the same scaling exponent, 

i.e., ip(t) oc but the probabilities of the leftward and rightward walk are different. 

Reduced map. — As can be seen in Fig. 1, points near the fixed points, x = and x = 1, 

on [0, 1] move to the left neighboring cell or the right neighboring cell, respectively. By 

translational symmetry, we can reduce an orbit of the map T(x) to that of an intermittent 

map on [0, 1] (Fig. 1). For example, we can obtain the following reduced map, 

X -|- (^) , mod 1 X < c, 
R{x) = I (4) 

X — {jz^y , mod 1 X > c. 



with < c < 1. The invariant density of the reduced map is given by 



p{x) = h{x)x 
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(5) 



where h{x) is a continuous function satisfying h{0) ^ and /i(l) 7^ [25|. Thus, the 
invariant density cannot be normahzed for 2; > 2. Consider the observation function 



f{A = { 



-1, X e [o,ci), 

0, X e [ci,C2) 

+1, X e [c2, 1), 



(6) 



where -R(ci) = 1 (ci < c) and i?(c2) = (c2 > c). It follows that X„ = /(xi) + ■ • • + /(x^) 
is regarded as a one-dimensional random walk, where x„, = R^{xq). The random walk X„ 
corresponds to a Levy walk where the persistent times distribution obeys a power law with 
exponent (3 = z/{z — 1). 

Dependence of EAMSD on ensemble. — Ensemble-averaged MSD (EAMSD) and TAMSD 
are defined by 
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{xl)= hm 1 ^(T™^-'^^ 

K^oo K ^ 
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where x^ is the kth initial point, and 



Sxl, 
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n— >-oo 77, 



(8) 



fc=0 



where fm{x) = {T"^{x) — x)^, respectively. For a finite invariant measure (2; < 2), ergodicity 
holds, i.e., EAMSD = TAMSD if initial points are distributed according to the invariant 
density of the reduced map. We note that ergodicity does not hold when the invariant 
measure cannot be normalized (infinite invariant measure) because an equilibrium ensemble 
cannot be reproduced. The impossibility to reproduce an equilibrium ensemble leads to 
aging 26|. In unbiased cases (c = 0.5), EAMSD is studied by the renewal theory and 



continuous time random walk l23l \2 



{xI^)e oc < 



^3-1/(2-1)^ 3 

mlnm. 



^ < z <2, 
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(9) 
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For c < 0.5, drifting motion arises in a direction toward the right. 



Unhke hyperbohc maps, statistical quantities determined by ensemble averages signifi- 
cantly depend on an initial ensemble in intermittent maps. In particular, it is shown that 
the behavior of the correlation function and the power spectrum density depend on an ini- 
|. In renewal theory 29(|, there are two well-known processes, i.e., ordinary 



tial ensemble 



renewal and equilibrium renewal process. An initial ensemble corresponding to a specific 
renewal process is reproducible in dynamical systems. In particular, an initial ensemble 
for an equilibrium renewal process is an absolutely continuous invariant measure. Figure 2 
shows that the EAMSD depends on an initial ensemble. If an initial ensemble of EAMSD 
has an invariant density (as in an equilibrium ensemble), all TAMSDs are equal to EAMSD 
(see Fig. 2). 



Lamperti- Thaler's generalized arcsine law. — In dichotomous stochastic processes an, an 
observable determined by the time average of an observation function g{an) is known to show 
random behavior if the mean residence time of a state diverges. In particular, the ratio of 
the occupation time of a state, Nn/n, does not converge to a constant, but converges in 
distribution, where A^„ is the occupation time of a state up to time n. The most classical 
example is the arc-sine law in coin-tossing: the distribution of the ratio of the period that 
a player is on the positive side converges to the arc-sine distribution 30|]. Lamperti showed 



that the distribution of the ratio of the occupation time of a state converges to the generalized 



arc-sine distribution in general dichotomous stochastic processes 
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In a dynamical system, the divergence of the mean residence time implies that the in- 



variant density cannot be normalized 



32| . Thaler has shown that the distribution of the 



;ime average of a characteristic function converges to the generalized arc-sine distribution 
8|-[l0|]. The generalized arc-sine law is valid for the time average of a non-L^(u) func- 
tion, i.e., J \f\dn = oo ll|. Lamperti-Thaler's generalized arc-sine (LTGA) law |9| states 
that for a map S : [0,1] ^ [0,1] satisfying (i) S{[0,c]) = [0,1] and S{[c,l]) = [0,1], (ii) 



S'{x) > 1 on (0,c] U [c,l); S'{0) = S'{1) = 1, and (iii) S{x) - x 



aox 



X — )■ 0) and 



X 



S{x) 



X 



X — 1) with p > 1 and constants aQ,ai > 0, and c G (0, 1), the 



time average of the observation function g{x) with g{0) = a and g{l) = b converges in 



distribution: 



Pr 



k=0 



where a = 1/p, 
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^ 1 - (^) (a < b) 
S'{c^) 



(10) 



{ao/a,y/PS'{c^y 
and the probabihty density function (PDF) is given by 

/3sina f^-^l-t^-^ 
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(11) 



(12) 



/32t2« + 2/3t"(l - t)" cos vra + (1 - t)^ 

This distribution is called the Generalized arc-sine distribution, which emereres in a subd- 
iffusive t..a„spo.t fl a.d weaU.y „o,.e.,od,c sta.is.ica. pKys.cs Q. The .ea„ of the 

occupation time for a state < c, (Nn/n), is given by (Nn/n) = 1/(1 + /3). The ensemble 
average of the time average of g{x) is given by (a + 6/3) /(I + /3). We note that the exponents 
a and /3 are determined by a behavior near the indifferent fixed points. 

Distributional response to a bias. — LTGA law cannot be applied to the reduced map of the 
asymmetric deterministic diffusion model with z >2 straightforwardly because the reduced 
map does not satisfy the condition (i). However, the condition (i) is not crucial because 
an important point in LTGA law is the reinjection to the indifferent fixed points. In fact, 
the conditions (i), (ii), and (iii) in {q] was generalized in 10|. The reinjection to the fixed 
point X = is determined by \imx->ci+o R'{x) and limx_>.c2+o R'{x). Moreover, the reinjection 
to the fixed point x = 1 is determined by lim^-^ci-o -^'(2^) and lim2;^c2-o -^'(2^)- Because 
of lima.^ci+o-R'(a;) = \ircix^ci~o R' (x) , lim^-^ca+o -R'(a;) = lim.x^c2-o R' i^) and R'{x) > 1 on 
(0, Ci] U [ci,C2] U [c2, 1), the way of a reinjection to x = and x = 1 is the same. Since 
the behavior near the indifferent fixed points is given by ^(x) — x ~ (x/c)^(x — )■ 0) and 
X — S{x) ~ {(1 — x)/(l — c)}^(x — )■ 1), we can apply LTGA law to the reduced map i?(x). 
Then, the exponents a and /3 are given by 
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z-1 



and /3 



(13) 



First, we consider ' 
tion with finite mean 



AMSD, where the observation function fm{,x) is the Lj^cfii^^ 1) func- 
ll|. By LTGA law, /„(0) = and /„(1) = m^, TAMSD converges 



to m^: 6x1, = for all c and z > 2. We note that the ballistic behavior, 6x1, = m^, is not 



due to a drift for c = 0.5. Next, we consider the time-averaged drift (TAD) defined by the 
time average of Vm{x) = T^{x) — x: 

n-l 



5xm = hm 



n— >oo n 

k=0 



^^t;^(xfc), (14) 



which is also the Lj^cfii^'' ^) function with finite mean. Figure 3 shows TADs obtained from 
different trajectories. By LTGA law, Vm{0) = —m and VmiX) = rn, we have 



Pr {5x„,/m <t}^l- (^) 



(15) 



This distributional limit theorem states that the time-averaged velocity, defined by V 



Sxm/'m, under a bias is intrinsically random, i.e., distributional response. Numerical simu- 
lations are in good agreement with the theory (see Fig. 4). The ensemble average of V is 
given by 

V^(F). = i^. (16) 

where {■)f is the ensemble average under an external bias. Figure 4 shows the response of 
to a bias c. 

Generalized Einstein relation. — The asymmetric parameter c with < c < 1 is con- 
sidered to be the probability of leftward walk if a persistent motion is terminated. Let 
c = e~'^ I {e~^ -|- e*&) be the leftward walk probability, where F is an external force, T is a 



temperature, and k is the Boltzmann constant 



22|. We consider a small bias, c = 1 /2-F/kT 



for F — > 0. Expanding /3 around c = 1/2 and substituting it into V, we obtain the linear 
response of V to a bias F: 

Although TADs are intrinsically random, the TAMSDs are not random and grow as for 
z > 2. Therefore, we have the following generalized Einstein relation for super diffusion: 

z F 



{6xr^)^ = -2—-^ — ^6xl. (18) 

Discussion. — We have found a distributional response in deterministic superdiffusion us- 
ing the distributional limit theorem in infinite ergodic theory. In a recent study we obtained a 
generalized Einstein relation for single trajectories using Hopf's ergodic theorem jssj. How- 
ever, this theorem does not work in deterministic superdiffusion because the observation 
function Vm{x) is not an L^{fJ,) function. Therefore, a distributional response is essential in 



super diffusion. Moreover, we noted that the generahzed Einstein relation ( ITSjl is different 
from that in anomalous diffusion We hope that our finding, i.e., the distributional 

response, will be observed in experiments of random time-averaged observables. 
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FIG. 2. (Color online) Mean square displacements {z = 1.9 and c = 0.5). Different symbols 
are the TAMSDs calculated by different initial points. Dashed and dotted lines are the EAMSDs 
based on equilibrium and ordinary ensembles, respectively, where an ordinary ensemble is a uniform 
ensemble on [0, 1] and an equilibrium ensemble is the points after 10^ times iterations. The slope 
of the solid line is the theoretical exponent ([9|). 
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FIG. 3. (Color online) Time averag ed drift [z = 3.0, c = 0.45, and n = 10^). TADs are calculated 
by different initial points. 
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FIG. 4. (Color online) Probability density function oiV{z = 2.5 and n = 10^). Numerical results 
are represented by the solid and dashed histograms with the theoretical PDFs (jlSp for c = 0.5 and 
c = 0.4, respectively. 
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FIG. 5. (Color online) Response to biases. Different curves are the response curves (|16p for 
different z. Symbols are the results of numerical simulations in the inset figure. 
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